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A thin layer of overdense plasma is created when an electromagnetic pulse interacts with a rapidly 
ionizing thin foil. This layer will reflect the incoming pulse, forming a so-called plasma mirror. A 
simple realistic model based on paired kinetic and wave equations is used to describe analytically the 
process of mirror formation and the reflection and transmission of the incident pulse. The model 
incorporates the exact description of the ionization process in the foil and the polarization and 
conduction currents that follow. The analytical description of the reflected and transmitted pulses 
as well as their dependence on foil parameters, and initial pulse amplitude and form are presented. 
Possible application and effectiveness of this process to improve laser pulse contrast are discussed. 
In the case of the linearly polarized incident pulse, there harmonic generation occurs, that is absent 
in the case of the circular polarization. The spectra of the reflected pulses for different initial forms 
and amplitudes are studied. 
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I. INTRODUCTION 

The interaction of ultrashort, high- intensity laser pulses with matter has lead to a new regime of nonlinear optics 
where the optical properties of the medium undergo ultrafast changes due to ionization produced by a laser pulse. 
In this context, by ultrashort, high-intensity pulses we refer to pulses with typical duration in the sub-picosecond 
regime, and whose intensity is high enough to cause field ionization of matter within a very few optical cycles, or even 
a sub-cycle interval. Therefore, by ultrafast changes of, e.g., a change in the index of refraction we mean that a latter 
varies significantly over a laser cycle, due to the rapid increase of the free electron density due to field ionization. As 
a consequence the pulse properties, such as envelope shape, spectrum and duration, can be strongly modified. 

A number of papers have addressed the problem of modeling the pulse propagation in a fast ionizing medium (see 
6-g- and studied the possible applications of ionization effects for pulse shaping and manipulation. The case of 
a high-density material, such as a solid target, is of particular interest since the free electron density produced by 
ionization can become quickly high enough to make the plasma refractive index turn from real to imaginary values, 
leading to self-reflection 0, y|. This effect is the physical basis for the so-called "plasma mirror" where plasma 
production at the surface of a mirror is used experimentally to obtain "clean" ultrashort pulses by removing the 
leading edge of the low- intensity "prepulse" or "pedestal" , which characterizes femtosecond laser systems d, H, H, 01 , 
by letting it be trasmitted through the mirror before the ionization-induced transition to opacity occurs. The resulting 
improvement of the laser pulse contrast (defined as the ratio between the prepulse energy and the short pulse energy) 
is important for various applications of current importance. These latter include, for example, the acceleration of 
fast ions from solid targets (see Q and references therein) where avoiding plasma formation at the target surface by 
the prepulse may enhance the ion source efficiency @. Similarly, high-contrast pulse may optimize high harmonic 
generation from solid target surfaces Another proposed application is the use of the combined effects of ultrafast 
optical shuttering and harmonic generation to obtain a sub-femtosecond pulse in transmission (llj . 

In these application regimes a thin (i.e., having a thickness less than an optical wavelength) foil target may be used 
in principle because the thickness of the reflecting plasma layer is usually of the order of the plasma skin depth. On 
the other hand the foil thickness may be used as an additional control parameter for applications. From the point of 
view of theoretical modeling a thin foil target is of interest because it allows an analytical treatment to some larger 
extent (see e.g. (Til [l^). thus providing an useful reference case for models and numerical simulations, and possibly 
aiding further optimization of the plasma mirror functionality. As a further reason of general interest, the laser- 
plasma interaction dynamics in a thin foil includes peculiar features that are not encountered either in underdense or 
overdense plasmas [ij, [l^ • 

In the present paper we utilize a simple analytic model for a dielectric foil ionization by an electromagnetic pulse. 
In order to describe the circularly polarized pulse interaction with the foil we use a kinetic equation paired with a 
wave equation. The effects of ionization are incorporated into the source term in the kinetic equation in the form 
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of a differential ionization rate [15|, |16[ . The solution of the kinetic equation, the momentum distribution function, 
enters the expression for the conduction and polarization currents that emerge in the foil. We point out differing 
behaviors between conduction and polarization electric currents. While the first grows with the increase of charged 
particle number, the second is proportional to the ionization rate, which means that it will be most important in the 
beginning of the ionization. In order to describe the field evolution during the interaction with the foil wc utilize 
the wave equation with nonlinear boundary conditions, obtained in Ref. [12l |. We show that as the pulse ionizes the 
foil its reflectivity grows giving rise to the reflected pulse. The form of the reflected and transmitted pulses depends 
strongly on the initial amplitude of the pulse. Wc study the process of plasma mirror formation under the action of a 
pulse with a gaussian envelope that has a pedestal. There is a significant pedestal reduction after the pulse reflection 
from the ionized foil. However there is a dependence of the reduction efficiency on the initial contrast of the pulse, 
which we study. We also consider the interaction of a linearly polarized pulse with a thin dielectric foil. In this case 
high order harmonics are generated due to the nonlinear dependence of the ionization rate on the pulse amplitude. 
Wc study the high order harmonic generation efficiency dependence on the initial amplitude of the pulse and on the 
form of the pulse. 

The paper is organized as follows. In section 2 we review the solution of the Maxwell equations with nonlinear 
boundary conditions. In section 3 we carry out the kinetic description of the currents that enter the wave equation. 
In section 4 we consider an analytical solution of the simplified wave equation. The results of the numerical solution 
in the case of a circularly polarized initial pulse are presented in section 5. In section 6 a linearly polarized pulse 
interaction with a foil is studied numerically. We summarize conclusions in section 7. 



II. WAVE EQUATION WITH NONLINEAR BOUNDARY CONDITIONS 

In this section we review the solution of the wave equation with nonlinear boundary conditions as it was presented in 
Rcf . [l^ . We consider the case of a circularly polarized plane wave impinging at normal incidence on a hydrogen-like 
foil, which is located in the plane a; = 0. The interaction of the pulse with the foil is taken to be uniform along the 
foil surface, which is justified by the fact that plasma mirrors operate in the regime in which the laser pulse is weakly 
focused on a mirror with a spot size of hundreds of wavelengths [1, H, 0, 0| • Relativistic effects of the ionization of 
the foil are neglected, which is valid if the magnetic field contribution to the ionization probability as well as to the 
electron motion is small. The density of the foil is taken to be high enough to ensure that the plasma density is not 
restrained by the fact that the number of neutral atoms might be exhausted. Also it is assumed that the electrons that 
emerge in the foil during the ionization remain localized within the foil due to the charge separation field produced 
in the presence of immobile heavy ions. So, the current is proportional to S{x), where 6{x) is Dirac delta function. 
The wave equation for the vector potential A{x,t) is 

duA - c^d^^A = 47r c (5(x)J(A) + &'{t)A{x, 0) + b(t)A{x, 0). (1) 

The first term on the RHS of Eq. ([1]) describes the electric current, localized in the foil, which is a functional of 
the vector potential taken at x = 0. The last two terms are equivalent to the initial conditions A(2;,0) = Ao(2;), 
(9tA(a;,0) = Ao(x), where Ao(a;) and Ao(a;) define the incident electromagnetic pulse. 
The solution of Eq. ([T|) can be presented as [l^ ■ 

t—\x\jc 

A(x-,i) = Ao(a;,t) + 27r ^ J(A(0,T))dr (2) 



By differentiating this equation with respect to time or spatial coordinate x wc obtain the explicit form of the electric 
or magnetic field respectively as functions of the incident electric field Eo(a;,i) and A(0,t), which can be found from 
Eq. © taken at a; = 0. 

E(x, t) = Eo(a;, i) - — J (A(0, i - |a;|/c)) , (3) 
c 

27r 

B(x, t) = Bo(a;, t) ^ e^. x J (A(0, t - |a;|/c)) sign(a;), (4) 

c 

where e^; is the unit vector along the direction of wave propagation. In order to solve these equations one should 
determine the expression for the current. 
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III. KINETIC DESCRIPTION OF THE CURRENT 



The behavior of the electrons in the foil in the presence of the homogeneous time-varying electric field can be 
described by the kinetic equation of the following form 



dfip,t) 
dt 



eE(i) 



df{p,t) 
dp 



g(p,E,t), 



(5) 



where the distribution function f{p,t) is normalized in such a way that / f{p,t)d^p gives the number of electrons 
per unit volume. (Similar approaches to field ionization modeling can be found, e.g., in Ref . p^l. [isj ) . Since we are 
considering the nonrclativistic regime, we can neglect the effects of the incident magnetic field, which are suppressed 
by the factor u/c, in the motion of electrons as well as in the ionization. The source term, q{p,'E,t), is connected 
with ionization of the foil by the incident electromagnetic pulse in such a way that J q{p, E, t)(fip gives the total 
ionization probability W ~ 8naUjKoK^£~^ exp [— 2k3/3£] of the unit volume in unit time (the time is measured in 
units of h? /me^). It is determined by the differential probability of ionization by a circularly polarized electric field 

M 



4k* UaUjKo 



■ exp 



2k3 
"3^ 



exp 



Pj^ £ (47ra)^c 
mc uj aB 



Kmc/ 



(6) 



where uj is the frequency of the incident electric field, £ = E/Ea is the normalized amplitude of the electric field, 
Ea ~ m?&' /h!^ = 5.14 • 10^ V/cm is characteristic atomic field, n = ^/TJIh, I is the ionization potential and 
Ih = me^ jif? ~ 13.6 eV is the ionization potential of the hydrogen atom, = I/huj, Ua is the density of 
neutral atoms, ub = hj^/me^ is Bohr radius, and a = 1/137 is fine-structure constant. The component of electron 
momentum perpendicular to the direction of pulse propagation is denoted by p±, and along the direction of pulse 
propagation is Px- The momentum distribution of produced electrons is very narrow with maximum at Px = 0, 
p±/mc = {£ /uj){16Tr'^ac/aB) p^ . [20| . This means that electrons preferably emerge from under the barrier with 
momentum perpendicular to the instantaneous direction of the electric field and equal to that of a free electron 
acquires in such a field. Since the momentum distribution is very narrow, we can use delta functions instead of 
exponents in order to simplify the calculations: 



■ exp 



2^3 
"3^ 



S[P^-^^^^]sipx). 



mc 



aB 



(7) 



We should also note that in the expression for the source term the absolute value of the field enters. 

We solve Eq.fO by integrating it along the particle characteristics. The equations for the characteristics are 



Po=Po(0), —— 
at 



(8) 



Introducing the function A.{t) 



J Eds, we obtain 





P± = -e Eds + p±o i-c, p± - eA{t) = p±o- 



(9) 



As a result the distribution function is given by the following functional 



z 

f^Mp.,P±-eA{t)]+ J q{p^-e[A{t)-A{t')],t'}dt', 



(10) 



where fo{px,P±) is the distribution fimction of the initial electrons before the passage of the electromagnetic wave. 
In our case /o — 0, since the foil was not ionized initially. The appearance of electrons due to the ionization of the 
foil gives rise to conduction and polarization electric currents. The conduction electric current is due to the electron 
motion in the field of the incident pulse, while the polarization electric current is due to the ionization process, since 
when the neutral atom is ionized an electric dipole is created. 
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In order to determine the form of the polarization electric current we write down the second moment of the kinetic 
equation 

(11) 

= E / [iKc^ + P^cY^^ - m^c^] q„(E, p, t)d^p, 



dK -s—^ dua f p_lc o 

^ -JcondE = S - 2^m„ — , icoud^e f—====dp, (12) 



which we rewrite as 

where K — "^J [(m^c^ +p^c^)^/^ — mcC^] fatPp is the kinetic energy and S = ^ /("^aC* + p^c^)^/^qa(E, p, i)d'^p. 

a Q 

Since in the foil not only electrons and ions are present but also neutral atoms, from which by the process of ionization 
the plasma is created, dn^/dt = dui/dt = —dua/dt. Here indexes e, i, and a stand for electron, ion, and atom 
respectively. Then 

2^"iaC = (me + m— ma)c — , (13) 

a 

where / — {nie + rrii — m^) c^, and is usually referred to as a mass defect or ionization potential, which is due to the 
fact that some electrons and ions are bound in atoms. In order to ionize an atom some amount of energy should be 
spent. So we rewrite Eq. (|12p in the following form 

^(K + In,) = j,o„dE + S. (14) 
at 

The energy balance equation for the circularly polarized electromagnetic wave is 

d /E^ \ 

( — ) = -jcoridE - jpo;E. (15) 



dt\AT: ) 

Here we have taken into account that no external electric current is present in our case. Adding these two equations 
we obtain the energy balance equation for the system of electromagnetic wave and particles 

d ( E^ \ 

A' + In, = -jpo,E + E. (16) 



The RHS of this equation should be equal to zero, because of the energy conservation law. Since we are considering 
a nonrelativistic case, we can neglect with respect to TOq. Then the expression for polarization current will take 
the following form 

_ E(t) dn, 

|E(i)|2' dt ■ ^ ' 

To our knowledge the polarization current was first introduced in Ref. j2l| while an extended discussion is given in 
Ref. [l3|. The expression given in Ref. [l3| is found to be identical to our formula (fT7|) when one realizes that the 
"ejection energy" of freed electrons £/ has been included in /. For our aims one may assume £/ = 0, i.e. we assume 
free electrons appear with zero energy, as an acceptable simplifying assumption; in such a case / represents in practice 
the standard ionization potential. In the more general case / becomes a smooth function of the electric field. 

It can be seen from Eq. (|17|) that the polarization electric current is along the direction of the electric field and 
is proportional to the ionization potential and to the ionization rate. Due to this fact we expect the contribution of 
the polarization electric current to be rather small. However, if we consider the ionization of heavy multicharged ions 
(where the binding energy of low lying levels is comparable to the electron rest mass) by an intense electromagnetic 
pulse, the polarization electric current can play an important role. It is similar to the case of electron-positron pair 
production by intense electromagnetic wave in plasma [2^ . where the polarization electric current is proportional 
to 2171,. In that case the polarization electric current totally changes the behavior of the wave, forcing the vector 
potential to oscillate around some nonzero value. Also the components of the electric field begin to oscillate with 
different frequencies. So in the process where the backreaction of the produced particles on the properties of the wave 
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that created them takes place, the polarization electric current should be taken into account. In the case of a very 
sharp electromagnetic pulse front, consider, for example, the interaction of high contrast and very intense pulse with 
the foil, the ionization rate is very high, and the polarization electric current will be significant, since the conduction 
electric current will not be able to develop large value due to the absence of electrons and ions in the first moments 
of evolution. 

Using the explicit expression for the distribution fmiction pop and the source term ([7|), and substituting ()12I17|) 
into Eq. we obtain: 



E(a::,t) =Eo(a:,i) 



2tt e^l 



\x\/c 



A{0,t-\x\/c) 



■ exp 



2k3 
"3^ 



dt' 



2n E(0,i- |a;|/c) SIln^KonaU 
'~ \E{0,t-\x\/c)\^ '£~ 



■ exp 



2^3"! 



(18) 



We can rewrite the equation for the electric field in terms of dimcnsionless field and vector potential, r/ = eE/mcwg, 
a = eA/mc^, where ujo is the frequency of the incident wave. The density of the plasma, that emerges due to the 
ionization, is taken in the units of hq = muj^^/ Aire^ ^ the critical plasma density, corresponding to the frequency 
of the initial pulse with time in the units of inverse frequency, ldq. The equation that governs the behavior of the 
dimcnsionless electric field is 



7?(a;,t) = riQ{x,t) - epa(0, i - \x/c\)n{Q,t~ \x/c\) 



I r;(0,t-|x/c|) 
'mc2 |r7(0,i- \x/c\)\' 



.n(0,i- |x/c|). 



where 



(19) 



n{t) = ^ 
no 



■ exp 



3\ri{t')\ 



dt' and n{t) 



no |?/(t)| 



■ exp 



(20) 



The normalized atomic field is denoted as rja ~ eEa/muJc and the normalized density of the foil, — uj'^Jj/2ll)qc, 
is a crucial parameter of the laser pulse interaction with a thin plasma layer introduced in Ref. [13[. It plays an 
important role not only in the process of electromagnetic pulse reflection/transmission at the thin foil [T^ . but also 
in the process of ion acceleration [l^ , when an intense laser pulse interacts with the thin foil, as well as in the process 
of high order harmonic generation in the relativistic regime of laser- foil interaction [23j . 

Equation ()19|) demonstrates that the conduction electric current is proportional to the number of electrons (ions), 
i.e. to plasma density, while the polarization electric current is proportional to the ionization rate, which makes it 
especially important in the first moments of intense ionization, should it occur. This distinction is due to the fact 
that the conduction electric current is connected with the motion of charged particles, while the polarization electric 
current is due to dipole creation through the ionization of neutral atoms. 

Due to the nonlinear dependence of the current on the electric field amplitude the above equation is difficult to 
solve analytically. 



IV. MODEL CASE 



In order to study the properties of Eq. ([19]) we consider a simplified model case. Let us assume that n{t) and thus 
n(t) are given functions of time. Then the equation p^. governing the evolution of the electric field, can be rewritten 
at the foil {x = 0) in the following form 

77 = ?7o-iVa+^iVi, where N = epuit) and A^i = e^-LM^. (21) 

The initial electromagnetic wave is 

VQ=mne'K (22) 
We shall look for the solution of Eq. ([2T|l in the following form 

rj = r]e'^^+^l (23) 
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After substituting the initial field and the solution into Eq. pT|) we obtain a system of algebraic equations for the 
amplitude and the phase of the transmitted wave: 



(24) 



(25) 



{rf' + A^i) cos ip = 7p]in — rj^N sin Lp 

{if + Ni)smip ~ ifN cos (p. 
Using the second equation we get the following expression for ip: 

rfN 

p = arctan — . 

Substituting this expression into the first equation of the system (HH), we find a forth order algebraic equation for 77: 

T^\l + N^) + rf{2Ni ~ rj^n) + = 0. (26) 

The solutions of this equation are 

1/2 

(27) 



r? = ± 



^(2iVi - r/lJ + V(2iVi - rjtJ^ - 4(1 + N^)N't 



2(l + iV2) 



Two other roots of forth order equation were omitted since 77 should be real. 

Let us consider the moment of time, when there is no longer any ionization, because all the atoms that could be 
ionized by the pulse have been ionized. The plasma density reached its maximum and the main part of the incoming 
radiation is reflected by plasma mirror. In this case Ni has dn{t)/dt = 0, and 



?7 



±- 



(28) 



This result states that even for plasma density much greater than critical there is a transmission of the incoming pulse 
through the ionized foil. However this transmission is suppressed by increased plasma density. 
Let us investigate the phase behavior. Using Eq. p5|) . we find 



The phase variation with time is 



p ~ arctan ■ 



dp 
'dt 



rfN 



N 

1 +iV2 



arctan N. 



> 0. 



(29) 



(30) 



So there is a blue shift in transmitted pulse frequency, while the density of plasma increases (see also [2J]). However 
in the region with N = const and Ni = there is no frequency shift and the radiation is transmitted without any 
frequency change. This fact should lead to the steepening of the transmitted pulse front and the rear should be of the 
same (half)width as the initial pulse. This can be seen from the results of numerical calculations, presented in the 
next section of this paper. So here unlike the case of the wave propagation in the ionized medium, the blue shift is in 
effect only when the ionization occurs. The frequency up-shift is connected with the rate of plasma density increase. 

V. THE RESULTS OF NUMERICAL SOLUTION 

In this section we present the results of numerical solution of Eq. at the foil, x = 0. We consider two cases of 
initial pulse form. First, the case of the pulse with gaussian envelope and gaussian pedestal, the maxima of the pulse 
and the pedestal temporaly coincide, 



E{0,t)=Eo{pcxp 



t-T/2 

^72~ 



+ (1 - p) exp 



t-T/2 

r/20 



iujQt 



(31) 



where p is the ratio of pedestal amplitude to the main pulse amplitude and r = 40 T is the pulse duration, T = 2t: /uj^). 
Second, the case of 



£:(0, t) = £'o[sin(7rt/r)]^ cxp[iwo^] 



(32) 
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a) b) c) 

FIG. 1: Evolution of the reflected (dotted curve) and transmitted pulse electric field (a), plasma density (b) and polarization 
current (c) for initial pulse of the form (|3H) . Time is measured in units of wave period {2n/uj), field is measured in units of 
mcuj/e, density is measured in units of critical density, Ucr, and currents are dimensionless. The parameter SpKo — Ax 10'^, 
X = 0. The first row corresponds to the initial pulse amplitude equal to rj^"'^ = 0.05, the second row to lyj""^ = 0.2. 



We consider time interval (0,t) to ensure that only one pulse is taken into account. This form of the pulse allows us 
to investigate the fronts of the reflected and transmitted pulses, since at i = the initial field is equal to zero. 

In Figs. 1 and 2 we present the results of numerical solution of Eq. (fT9|) with initial pulse having form ((3T|) and ((32|) 
respectively. In Fig. 1 maximal electric field is 77™"^ = 0.05 (first row) and 77™°^ — 0.2 (second row). The amplitudes 
of the reflected and transmitted pulses arc shown in Fig. la. The plasma density evolution which is closely connected 
with the conduction current is shown in Fig. lb. The behavior of polarization current is presented in Fig. Ic. In 
the case of ry™"^ = 0.05 (Fig. 1, first row) the incoming pulse is not intense enough to sufhcicntly ionize the foil, the 
resulting plasma density is 0.5 of critical. Due to this fact the reflectivity of the of the ionized foil is low, and the 
reflected pulse is much lower than the transmitted one. However it can clearly be seen that the pedestal is transmitted 
through the foil. The reflected pulse starts to differ from zero only when the main part of the initial pulse arrives at 
the foil. 

In the case of 77™"^ = 0.2 (Fig. 1, second row) the situation is different. Since the initial pulse intensity is high 
enough to produce a plasma with a density 4 times larger than critical by means of foil ionization, the reflectivity of 
the foil tends to unity, as shown in Fig. la the main pulse is almost completely reflected. However we should note 
that the intensity of the pedestal is not enough to ionize the foil, that is why it is completely transmitted through 
the foil. When the main pulse arrives at the foil it quickly ionizes the foil, it can be clearly seen in Fig lb, second 
row, where the plasma density evolution is presented. In this case the reflected pulse is pedestal-free. So through the 
interaction of the intense electromagnetic pulse with rapidly ionized foil the reflected pulse is almost the same as the 
initial one apart from the fact that the pedestal is greatly reduced in the reflected pulse. We should also note that 
the reflected pulse has a steeper front than the initial one and its tail is of the same (half)width as the initial one (see 
Fig. la). These features are qualitatively expected from the analysis of the simplified wave equation carried out in 
the previous section. 

In order to illustrate the contribution of the polarization current to the evolution of the reflected and transmitted 
pulses, we present Fig. Ic. As it was pointed out above, the expected polarization current contribution should be 
rather small, since it is proportional to the ionization potential. However the polarization current is also proportional 
to the ionization rate, so it can be important in the first moments of intense ionization. It can be seen from Fig. 
Ic that the polarization current is small compared to conduction one (to estimate the amplitude of the conduction 
current; one can multiply electric field amplitude by the number of particles), but it is non-zero in the region where 
the conduction current is far from its maximal value, and thus its contribution can be of the order of conduction 
current contribution. 

In Fig. 2, where the results of numerical solution of Eq. for initial pulse (|32p are presented, maximal electric 
field is rj^'^^ = 0.05 (first row) and 77™°^ = 0.5 (second row). The amphtudes of the refiected and transmitted pulses 
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FIG. 2: Evolution of the reflected (dotted curve) and transmitted pulse electric field (a), plasma density (b) and polarization 
current (c) for initial pulse of the form (|32p . Time is measured in units of wave period (27r/a;), field is measured in units of 
mcuj/e, density is measured in units of critical density, Ucr, and currents are dimensionless. The parameter €pKo — Ax 10'^, 
X = O.The first row corresponds to an initial pulse amplitude equal to r/JJ^"'^ = 0.05, the second row to 77™"^ = 0.5. 

are shown in Fig. 2a. The evohition of plasma density can be seen in Fig. 2b, and the behavior of polarization 
current in Fig. 2c. In the case of the initial pulse maximal amplitude of 77™"^ = 0.05 (first row), the reflected and 
transmitted pulses are roughly of the same order, but the reflected pulse has a steeper front, since the reflectivity of 
the foil depends on the density of plasma, which evolution is shown in Fig. 2b. With increasing of the initial pulse 
amplitude the situation changes, in the case of 77™°^ = 0.5 (second row) the transmitted pulse is much smaller than 
the reflected, which almost coincides with the initial one. Since the amplitude of the initial pulse is high enough to 
create a significant amount of plasma even in the first moments of the pulse interaction with the foil, the initial pulse 
gets reflected. Here, analogous to the case of the pulse with the pedestal in interaction with the foil, the reflected 
pulse has the characteristic temporal profile. Its front is steeper than its tail, while the latter (half)width is the same 
as the (half )width of the initial pulse. The properties of the polarization current are the same as in the previous case 
of pulse with pedestal. The absolute value of the polarization current is much smaller than that of the conduction 
current. However the polarization current has a maximum in the region where the conduction current is far from 
reaching its maximal value. 

Let us consider the case of the pulse with pedestal once more but from the point of view of obtaining high contrast 
pulses. In order to get a high contrast pulse the pedestal should be transmitted through the foil and the safe-reflection 
should be triggered by the main pulse. In Fig. 3 we present a contour plot of the reflected pulse resulted contrast 
dependence on the initial pulse contrast and intensity. It is calculated as the ratio of the reflected pulse pedestal 
instantaneous intensity at the point where the pedestal and the main pulse envelopes intersect and the reflected pulse 
maximum intensity. The contours correspond to the values 10^*"', k = 1..10 of the resulted contrast. We see that for 
pulses with the low initial contrast there is no way to improve it by the reflection on the rapidly ionized foil. It is 
due to the fact that the intensity of the pedestal is high enough to ionize the foil giving rise to the reflection of the 
incoming radiation. When the initial contrast improves there appears a region (to the right of the dashed line in Fig. 
3) of initial pulse intensity values where the resulting contrast of the reflected pulse is better than that of the initial 
pulse. With the further improvement of initial pulse contrast this region enlarges, and for initial contrast better than 
10~^, this region covers all the range of considered amplitudes of the initial pulse. So in order to improve the laser 
pulse contrast in the interaction with ionizing foil the initial pulse contrast should be high enough not to allow the 
prepulse trigger the foil ionization-induced transition to opacity. 
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FIG. 3: Dependence of the reflected pulse contrast (the ratio of the pedestal intensity to the main pulse intensity) on the initial 
pulse contrast (horizontal axis, logarithmic scale) and main pulse dimensionless amplitude (vertical axis, the fleld is measured 
in units of mcuj/e). The contours represent the values 10"^'', k — 1..10. The dashed curve depicts values where the resulting 
contrast is equal to the initial one. To the right of this curve lie the values of the improved contrast. 

VI. HIGH ORDER HARMONIC GENERATION BY THE LINEARLY POLARIZED PULSE. 

Since a circularly polarized pulse docs not generate high order harmonics when it ionizes the foil, we consider the 
ionization of a foil by a linearly polarized pulse. In this case all the formulae of sections 2 and 3 are valid, the only 
thing that should be taken into account is that in this case the electrons emerge from under their barriers with a 
momentum distribution that has its maximum at zero. However this distribution is very narrow [Toj . its width is 
much less than the typical momentum that an electron acquires in the electric field. Due to this fact the initial 
momentum distribution can be approximated by Dirac delta function in the expression for the source term. In this 
case the equation (fT9)) should be modified, and it becomes 



In what follows we present the results of numerical solution of equation (|33p for the same initial pulses as in the 
preceding section. In Figs. 4 and 5 the form of the reflected and transmitted pulses, plasma density evolution and 
the spectrum of the reflected pulses are shown for the initial pulse shapes ([31]) and ((32)) with the same amplitude and 
duration as in Figs. 1 and 2. The behavior of the reflected and transmitted pulses is almost the same as in the case 
of the circularly polarized initial pulse as can be seen from Figs. 4a and 5a. However the evolution of the plasma 
density is slightly different. In order to illustrate this we present Figs. 4b and 5b, where a step-like function is shown, 
its behavior is similar to that in the case of circular polarization apart from the step-like feature. This property of the 
density evolution is due to the ionization. Since the ionization rate is determined by the absolute value of the electric 
field, in the case of the linear polarization it has two maxima per period of the light wave. The step-like function just 
demonstrates that the ionization occurs when the electric fleld amplitude reaches its maximum, while there is almost 
no ionization in between. 

In Figs. 4c and 5c the spectra of the reflected pulses are presented. We see that due to the nonlinear dependence 
of the ionization rate on the absolute value of the fleld harmonic generation occurs. Since we consider only normal 
incidence of the laser pulse on the foil the electric field is parallel to the foil, and that is why only odd harmonics are 
generated as it can be seen from Figs. 4c and 5c. In order to compare the cases of linear and circular polarizations 
we also present in Figs. 4c and 5c the spectra of circularly polarized reflected pulses, that have the same energy as 
the linearly polarized ones. Here we can clearly see that there is no high order harmonic generation in the case of 
circularly polarized pulses. Whereas the modification of the spectrum that is due to the ionization blueshift, that was 
mentioned in Section IV, manifests itself in both cases. However this effect is weak when the plasma mirror works 
in the condition close to the optimal regime as the transparency-opacity transition is abrupt, since the interaction 
length is very small, of the order of the skin depth. 

Let us estimate the collective efficiency of high order harmonic generation. In order to do so we define the parameter 



t-\x\/c 




(33) 
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FIG. 4: Evolution of reflected (dotted curve) and transmitted electric field envelopes (a) and plasma density (b) in the case of 
linear polarization, the spectra of the reflected pulses in the case of linear (upper curve) and circular (lower curve) polarizations 
for the same pulse energy (c) for initial form of the pulse (|3ip . Time is measured in units of wave period (27r/ti^), field is 
measured in units of mccj/e and density is measured in units of critical density, Ucr- The parameter epKo = 4 x 10''. The first 
row corresponds to the initial pulse amplitude equal to 77™"^ — 0.05, the second row to r;™"^ = 0.1. 

R: 

UJ — OC 

J El^{u:)du: 

R-"^ ■ (34) 

This parameter has the meaning of the energy that is contained in higher harmonics in the units of initial pulse energy. 
In other words it measures the efficiency of initial pulse energy transformation into high frequency radiation. 

The dependence of R on the initial pulse amplitude is shown in Fig. 6 for different forms of the initial pulse: a) 
gaussian pulse with pedestal (j3ip and b) pulse with sine squared envelope (|32[) . The behavior of the curves is almost 
the same, apart from the fact that in Fig. 6a the parameter R is measured in the units of 10^'^, while in Fig. 6b in 
the units of 10~^. For small values of the initial pulse electric field amplitude R is equal to zero. This corresponds to 
the fact that there is no ionization for low values of Eq. Then at Eq « 0.1 there is a maximum of R. At this value 
of the initial field amplitude the reflected and transmitted pulses are roughly of the same magnitude. So the high 
harmonic generation due to foil ionization is most effective where the duration of ionization is the longest. 

VII. CONCLUSIONS 

A model for the interaction of a short laser pulse with a thin foil target undergoing field ionization has been 
derived. The ionization dynamics have been described starting from a kinetic equation for the electrons where the 
source term includes the ionization rate and the momentum distribution of freed electrons. Energy conservation has 
been included via a polarization current term. Using the electromagnetic boundary conditions at the foil surface a 
set of ordinary differential equations has been derived. Such a compact model can be used for simple analysis and 
preliminary parametric studies of the "plasma mirror" properties of the thin foil, such as prepulse suppression and 
ultrafast optical shuttering. As an example the transmission, reflection and shaping of a short laser pulse with a 
pedestal have been studied as a function of the initial pulse contrast and peak intensity. The generation of high 
harmonics by ionization has also been studied for a linearly polarized pulse. 
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FIG. 5: Evolution of reflected (dotted curve) and transmitted electric field envelopes (a) and plasma density (b) in the case of 
linear polarization, the spectra of the reflected pulses in the case of linear (upper curve) and circular (lower curve) polarizations 
for the same pulse energy (c) for the initial form of the pulse (|32p . Time is measured in units of wave period {2-n/uj), field is 
measured in units of mcui/e and density is measured in units of critical density, ricr- The parameter e^Ko = 4 x 10'^. The first 
row corresponds to the initial pulse amplitude equal to t)™"^ = 0.05, the second row to ry(™"^ = 0.5. 

7? X 10^ Ry. l(f 

5.0 p 




£^0 ' ' Eq ' 

a) b) 

FIG. 6: Effectiveness of harmonic generation vs the initial field amplitude for the cases of initial pulse form pip (a) and (|32p 
(b). The parameter epKo = 4 x 10^. 



Further extension of this compact model will require the inclusion of multiply ionized species and off-normal 
incidence laser pulses that can be achieved by the modification of the source term and the proper transformation of 
the reference frame correspondingly. This will allow the study of the "plasma mirror" operation with respect to the 
angle of incidence and different states of polarization. 
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